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TITLE:  NONLINEAR  AEROELASTICITY  OF  COMPOSITE  STRUCTURES 
(AFOSR  Contract  No.  F49620-86<C-0066) 

Principal  Investigator:  John  Dugundji 
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Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


OVERVIEW: 


There  are  few  analytic  studies  of  nonlinear  stall  flutter,  most  studies  being 
experimental,  and  recently,  of  the  computational  fluid  dynamic  type.  The  present 
investigation  develops  an  analytic  method  to  include  the  nonlinear  aerodynamic  and 
structural  effects  into  a  full  3-dimensional,  aeroelastic  problem,  using  the 
mathematical  tools  of  Fourier  analysis,  harmonic  balance,  and  the  Newton-Raphson 
method  as  a  numerical  solver.  The  nonlinear  aerodynamics  are  introduced  using  the 
ONERA  stall  flutter  model  ^  while  the  nonlinear  structures  are  introduced  by  moderate 
deflection  cantilever  beam  theories  such  as  Hodges  &  Dowell^  which  is  used  in 
helicopter  theory  and  involves  coupling  of  traditional  bending  and  torsion  modes 
with  fore-and-aft  modes.  It  is  hoped  that  this  new  nonlinear  analysis,  and  the 
attendant  flutter  experiments  on  various  aeroelastically  tailored  composite  wings, 
will  shed  important  light  on  the  phenomenon  of  high  angle-of-attack  stall  flutter,  its 
severity,  and  how  it  evolves  out  of  small  amplitude,  linear  flutter  at  low  angles-of- 
attack. 


SUMMARY  OF  ACCOMPLISHMENTS: 

A  combined  experimental  and  analytic  program  was  undertaken  to  investigate 
the  nonlinear  stall  flutter  and  divergence  behavior  of  aeroelastically  tailored 
composite  wings.  Experimentally,  a  set  of  zero  sweep  cantilever  wing  models,  of 
length  30  cm,  were  constructed  from  flat  composite  plates  covered  with  styrofoam 
fairings  to  give  NACA  0012  contours.  The  wings  had  varying  amounts  of  Ending- 
torsion  coupling  ranging  from  positive  to  negative  values.  Tests  were  conducted  in 
the  wind  tunnel  at  varying  root  angles  as  shown  in  Fig.l.  These  tests  revealed  linear 
flutter  and  divergence  as  well  as  nonlinear  torsion  stall  flutter  and  bending  stall 
flutter.  Flutter  velocity  boundaries,  frequencies,  limit  cycles,  and  static  deflection 
were  observed  at  different  root  angles  of  attack.  Analytically,  a  new  nonlinear 
flutter  analysis  technique  was  developed  based  on  the  ONERA  aerodynamic  stall 
model,  together  with  simple  frequency  and  harmonic  balance  techniques  to  obtain 
the  resulting  nonlinear  stall  limit  cycles  that  can  arise  from  aeroelastic  instabilities. 
Both  torsional  stall  and  bending  stall  flutter  limit  cycles  were  successfully  predicted 
by  this  technique  and  matched  the  observed  experimental  behavior.  Figure  2  shows 
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the  flutter  boundary  versus  root  angle  of  attack  for  a  typical  wing  with  no  bending  - 
twist  coupling,  while  Fig. 3  shows  the  boundary  when  positive  coupling  is  present. 
Accompanying  static  positions  and  limit  cycle  amplitudes  for  Fig. 2  are  shown  in 
Fig. 4,  indicating  a  torsional  stall  flutter  situation,  while  the  corresponding  quantities 
for  Fig. 3  indicated  a  bending  stall  flutter  situation  (also  note  the  much  lower  flutter 
frequency  in  Fig. 3  corresponding  to  the  bending  mode  rather  than  the  torsion 
mode).  A  detailed  report  on  this  work  is  given  by  Dunn^,  while  a  paper  summarizing 
some  of  these  results  was  presented  at  the  31st  AIAA  SDM  Conference  last  year^. 

An  additional  set  of  composite  wings,  twice  the  size  of  the  previous  wings  was 
also  constructed  and  tested.  This  yielded  improved,  high  quality  experimental  data  to 
compare  with  the  analysis,  and  to  observe  the  transitions  from  classical  linear  low 
angle  of  attack  flutter  to  stall  flutter  at  high  angles  of  attack.  Supplementing  this 
new  data,  an  expanded  nonlinear  analytical  effort  was  pursued  using  more  modes  and 
including  fore-and-aft  motion  of  the  wings  as  well  as  the  usual  bending  and  torsion 
motions.  These  fore-and-aft  motions  couple  in  nonlinearly  into  the  analysis,  and 
were  also  observed  during  the  tests. 

To  summarize,  an  analytical  method  was  developed  to  incorporate  nonlinear 
structural  and  aerodynamic  stall  effects  into  a  full,  3D  aeroelastic  flutter  analysis. 
High  quality  stall  flutter  data  was  obtained  to  check  the  new  nonlinear  analysis. 
Transitions  from  linear,  coupled-mode,  bending-torsion  flutter  to  single-mode, 
torsional  stall  flutter  was  observed  experimentally  and  predicted  analytically.  Also 
transition  from  linear  divergence  to  bending  stall  flutter  was  observed 
experimentally  and  predicted  analytically.  These  experimental  and  nonlinear 
analytic  results  should  be  of  interest  in  understanding  stall  flutter  at  high  angles  of 
attack. 
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FOREWORD 


The  following  pages  represent  a  more  detailed  summary  of  the 
work  performed  on  this  contract.  They  are  from  a  technical 
presentation  given  on  October  10,  1991  at  the  Air  Force  Office  of 
Scientific  Research  Contractor's  Meeting  on  Structural  Dynamics  at 
Datyon,  Ohio. 
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OBJECTIVES 


Investigate  nonlinear  structures  and  nonlinear 
aerodynamics  in  large-amplitude,  high  angle-of-attack, 
stall  flutter  of  aeroelastically  tailored  wings 

Develop  method  to  incorporate  nonlinear  effects 
(aerodynamic  and  structural)  into  a  flutter  analysis 


Develop  In  context  of  simple,  modal  JluJ^  4na^ 
to  keep  computational  costs  low 

Concurrently  develop  experimental  base^^^of  small  and 
large  amplitude  flutter  data  for  variety  of  composite 
laminate  wings 


ANALYSIS  -  RAYLEIGH-RITZ  MODEL 


w(x,y,t)  = 


(t) 


^22^1,  yyYj^yy  +  ^’^66Yi,xyYj,xy 

°12  ^Yi^x.-Yj,yy  Yi^yyYj^xx^  ^Yi^xxYj,xy  Yi^xyYj,xx^  ^ 

^^26  ^"^i/yyYj,  xy  '*'  Yi,xyYj,yy^  ^  ^ 


■JJ 


mY^Yj  dA 


Qi  =  Jjf  (x,y)  Y.  (x,y)  dA  =  J|f  (x,  y)  (j).  (x)  \|/^(y)  dA 


[M] {q}  +  [K] {q}  =  {Q} 

-  Mass  and  stiffness  of  styrofoam  calculated 
and  added  to  overall  matrices 

~  Empirical  cubic  stiffening  added  to  torsional  mode 

K22  =  K^2  "  ^^^2 
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ANALYSIS 


AERODYNAMICS 


if  "k  it 

C  z2  +  2dwCz2 


+  w2 (l+d2) Cz2  = 
=  -w2  (i+d2)[^ACz 


where, 


-  Aerodynamic  model  from  ONERA 

-  Coefficients  determined  semi-empirically 


-  2D  aerodynamics  at  each  spanwise  location, 
corrected  for  finite  span  (3D)  effects, 
integrated  to  give  modal  forces 


ANALYSIS 


FOURIER  ANALYSIS 


ANALYSIS 


FOURIER  ANALYSIS 


-  Assume  polynomial  form  to  aerodynamic  curve, 

Ji 

ACz(a)  =  ^a.j(a-ai)3  ;  ai<a<ai+i 

j=0 


-  Sinusoidal  input  for  effective  angle  of  attack, 

•  a(X)  =  CXo  +  avsin9  =  Oto  +  assin(kT)  +  ttccos  (kx) 
(P  =  kX  +  ^  =  -\/a2  +  a2  ;  ^  ^ 


-  Fourier  analysis  for  harmonic  components  of  AC; 


+7t  ^ 


=  -^JaCj  (X)sin9  d9 


EXPERIMENT 


Hercules  AS4/3501-6  graphite/epoxy  test  specimens 
[02/90]s,  [+152/0]s,  [-152/0]s  layups 


STYROFOAM  FAIRINGS 


EXPERIMENT 


Wind  tunnel  tests  in  5’  x  7.5’  free  jet  test  section 
Data  from  strain  gauges  and  visual  data  on  video 
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Flutter  boundary  and  frequency  variation 
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Flutter  boundary  and  frequency  variation 
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Tip  angle,  and  deflection  &  angle  oscillation  amplitudes 
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Flutter  boundary  and  frequency  variation 


RESULTS  -  [-152/0]s 


vaoarr  mm 


Tip  angle,  and  deflection  &  angle  oscillation  amplitudes 


Velocity  (m/s) 


[Og/QOJg  Flutter  boundary 


n  =3,  n  =3,  n  =3 

b  t  f 


5  10  15  20  25  30 

Root  AOA  (deg) 

Frequency  along  flutter  boundary 


Foam  E|^,  Ej,  C^j(MPa): 

20,  20,  10  based  on  static  deflection 
10,  10,  6  based  on  vibration 
15,  15,  8  compromise 
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Frequency  (Hz)  Velocity  (m/s) 


[  +  1  5^/ 0^]^  Flutter  boundary 


Frequency  along  flutter  boundary 


25-JUL-91 


9 


[-1  5  /  0  ]  Flutter  boundary 


to 


[0/90] 


I 


CONCLUSIONS 


Developed  analytic  method  to  incorporate 
nonlinear  structural  and  aerodynamic  effects  into 
a  full,  3D  aeroelastic  problem 

Method  used  with  simplified  aerodynamics  and 
lowest  harmonic,  but  can  be  extended  to  more 
complex  variations 

Experimental  data  obtained  to  check  analysis 


Transition  from  linear,  coupled-mode,  bending-torsion 
flutter  to  single-mode,  torsional  stall  flutter  observed 
experimentally  and  predicted  analytically 


Transition  from  linear  divergence  to  bending  stall 
flutter  observed  experimentally  and  predicted 
analytically 


